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ABSTRACT 

We present the results of the simultaneous deep XMM- Newton and Chandra 
observations of the bright Seyfert 1.9 galaxy MCG— 5-23-16, which is thought 
to have one of the best known examples of a relativistically broadened iron 
Ka line. We detected a narrow sporadic absorption line at 7.7 keV which 
appears to be variable on a time-scale of 20 ksec. If associated with Fexxvi 
Lyo; this absorption is indicative of a possible variable high ionization, high 
velocity outflow. The time averaged spectral analysis shows that the iron 
K-shell complex is best modeled with an unresolved narrow emission component 
(FWHM< 5000 km/s, EW ~ 60 eV) plus a broad component. This latter 
component has FWHM~ 44000km/s, an EW~50 eV and its profile is well 
described with an emission line originating from the accretion disk viewed 
with an inclination angle ~ 40° and with the emission arising from within a 
few tens of gravitational radii of the central black hole. The time-resolved 
spectral analysis of the XMM-Newton EPIC-pn spectrum shows that both 
the narrow and broad components of the Fe K emission line appear to be 
constant within the errors. The analysis of the XMM-Newton/RGS spectrum 
reveals that the soft X-ray emission of MCG-5-23-16 is likely dominated by 
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several emission lines superimposed on an unabsorbed scattered power-law 
continuum. The lack of strong Fe L shell emission together with the detection of 
a strong forbidden line in the 0 VII triplet supports a scenario where the soft X- 
ray emission lines are produced in a plasma photoionized by the nuclear emission. 


Subject headings: galaxies: active - galaxies: individual (MCG-5-23-16) - galax- 
ies: Seyfert - X-rays: galaxies 

1. Introduction 

One of the key issues in high energy research on Active Galactic Nuclei (AGN) is the 
study of the 6.4 keV iron Ka emission line profile, which can provide fundamental diagnos- 
tics of the physical and dynamical conditions of their central engines. The fluorescent Fe 
Ka emission line is a prominent and ubiquitous feature in the X-ray spectra of AGN and 
it is believed to be produced in the innermost regions of the AGN, such els the Broad Line 
Region (BLR), the circumnuclear obscuring torus and/or the accretion disk. The profile of 
the line itself provides primary information on the region from which it originates. If the Fe 
Ka emission line is produced far from the nucleus, e.g. in the putative torus, its profile is 
expected to be narrow, while if the line originates in the innermost part of the accretion flow 
a broad and asymmetric profile is predicted as a result of the special and general relativistic 
effects such as Doppler shifts, gravitational redshift and light bending (see Fabian et al. 2000; 
Reynolds & Nowak 2003 for a review). In the latter case the shape of the profile itself could 
be used to derive information on the nature of the black hole and accretion disk system. 
The observations with ASCA of relativistically broadened iron Fe Ka emission lines in AGN 
(Nandra et al. 1997) and in particular the detection of a broad and skewed profile in the long 
ASCA observation of MCG-06-30-15 (Tanaka et al. 1995) was considered the first evidence 
that at least some line emission originates in the inner part of the accretion disk close to 
the central black hole. However, the scenario emerging from XMM- Newton and Chandra 
observations of AGN appears to be more complex. Indeed, these observations have shown 
that only a handful of objects show the presence of the relativistically broadened line, while 
the narrow Fe emission lines at 6.4 keV is an ubiquitous feature in many type I AGN (see 
Bianchi et al. 2004; Reeves et al. 2004; Yaqoob & Padmanabhan 2004). Furthermore the 
broad component appears to be in general weaker than what was expected to be after the 
initial ASCA results and in some case it may be absent (i.e. NGC 4151,Schurch et al. 2003). 
These observations have also shown that the interpretation of the Fe profiles can be strongly 
dependent on the modeling of the underlying continuum, which can be complicated by the 
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presence of complex absorption (i.e. NGC 3783 Reeves et al. 2004; NGC 3516 Turner et 
al. 2005) and reflection components (see Reeves et al. 2007 and reference therein). Further- 
more, red- and blue-shifted Fe absorption lines, associated with the presence of infalling or 
outflowing matter in the proximity of the black-hole have been detected in the X-ray spectra 
of QSOs and Seyfert galaxies (see Cappi 2006 and references therein). These absorption 
and emission features together with the complexity of modeling the underlying continuum 
makes the study of the Fe line profiles more complex and thus feasible only for the brightest 
objects. 

In this framework MCG-5-23-16 represents one of the best and more robust examples of 
a relativistically broadened Fe line. MCG-5-23-16 is a nearby (z = 0.008486) Seyfert 1.9 
galaxy, with a typical 2-10 keV flux of ~ 8 x 10 -11 erg cm" 2 s _1 making it one of the X-ray 
brightest Seyfert galaxies. Previous X-ray observations showed that the X-ray emission of 
MCG-5-23-16 resembles the classical spectrum of a Compton thin (i.e. Nh < 10 24 cm _2 ) 
Seyfert 2, with a soft excess below 1 keV and a column density Nh ~ 10 22 cm -2 (Dewangan 
et al. 2003; Balestra et al. 2004). Higher energy observations (i.e. above 10 keV) performed 
with RXTE (Weaver et al. 1998; Mattson & Weaver 2004) and BeppoSAX (Perola et al. 
2002; Risaliti 2002a) detected a Compton reflection component, which was interpreted as 
reprocessed emission from the distant molecular torus. 

A strong broad Fe Ka line was first detected with ASCA (Weaver et al. 1997, 1998) with 
an equivalent width EW~ 200 eV, which could be modeled with a broad relativistic line 
profile ( i ~ 50°) plus a narrow core at 6.4 keV. The presence of both these components 
has been subsequently confirmed with Chandra and XMM -Newton observations. Two short 
XMM -Newton observations, whose summed exposure time was ~ 25 ksec, confirmed the 
presence of a underlying broad component with an EW~ 100 eV (Dewangan et al. 2003; 
Balestra et al. 2004). However the relatively short exposure time of the past XMM- Newton 
observations did not allow the above authors to put strong constraints on the geometry of 
the emission region. Indeed neither the inner disk radius nor the inclination angle could be 
accurately derived using only the XMM-Newton data. 

In this paper we present an analysis of the iron K line profile and variability from a simultane- 
ous deep XMM -Newton (130 ksec) and Chandra (50 ksec) observations of MCG-5-23-16; the 
analysis and results of the simultaneous deep Suzaku (~ 100 ksec) are described in Reeves et 
al. (2007). The XMM -Newton and Chandra observations and data reduction are described 
in section 2. In section 3 we present the modeling of the time-averaged XMM -Newton EPIC 
and RGS spectra and the results of the spectral fits of the simultaneous Chandra HETG 
spectra. In section 4 we report the results obtained with time-resolved spectral analysis, 
aimed to assess the possible variability of the iron K emission line during the long observa- 
tion and to investigate the relation (or lack of) between the Fe emission line intensity and the 
flux of the underlying X- r ay continuum. Finally in section 4 we discuss t he appea r ance of ~ a ~ 
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sporadic absorption feature at 7.7 keV (rest- frame) which is indicative of a possible variable 
high velocity (v~ 0.1c) outflow. The results are discussed and summarized in section 5. 


2. Observations and data reduction 

In December 2005, MCG-5-23-16 was observed simultaneously with many different X- 
ray observatories: Suzaku, XMM-Newton, Chandra and RXTE; in Table 1 we report the 
log of the different observations. In this paper we concentrate on the XMM-Newton and 
Chandra observations, while the Suzaku and RXTE observations and results are described 
in detail by Reeves et al. (2007). 


2.1. XMM -Newton 

MCG-5-23-16 was observed with XMM-Newton on 2005 December 8th for a total expo- 
sure time of about 130 ksec (see Table 1). The pn, MOS1 and MOS2 cameras had the medium 
filter applied; the MOS1 and MOS2 were operating in Small Window mode, while the pn 
was in Large Window mode. The XMM-Newton data have been processed and cleaned using 
the Science Analysis Software (SAS ver 6.5) and analyzed using standard software packages 
(FTOOLS ver. 6.1, XSPEC ver. 11.3). In order to define the threshold to filter for high- 
background time intervals we extracted the 10-12 keV light curves and filtered out the data 
when the light curve is 2 a above its mean. This screening yields net exposure times (which 
also includes a dead-time correction) of 96 ksec, 101 ksec and 103 ksec for the pn, MOS1 
and MOS2 respectively. For the scientific analysis we concentrated on the EPIC-pn data 
which have the highest signal-to-noise, and we used the MOSl and MOS2 data to check for 
consistency. Taking into account the brightness of the source (the 2-10 keV count rates are 
7.9 cts/s, 2.7 cts/s and 2.8 cts/s for the pn, MOSl, MOS2), we ran the sas task epatplot to 
check for possible pile-up and we found that both in the pn and MOS detectors the pile-up 
fraction is below 1%. However since we have good photon statistics, when analyzing the 
time averaged pn spectrum, we decided to use only the pattern zero data (which correspond 
to single events), which are better calibrated, and we considered the pattern 0-4 (single and 
doubles) when we extract spectra with lower exposure time for the time-variability study. 
The EPIC pn source spectrum was extracted using a circular region of 37" and background 
data were extracted using two circular regions with an identical radius (37") centered at 
~ 4' from the source. EPIC MOSl and MOS2 data were extracted using a source extraction 
region of 27" radius and two background regions with identical size (27") selected on the 
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nearby CCDs. Response matrices and ancillary response files at the source position have 
been created using the SAS tasks arfgen and rmfgen. Background subtracted data were 
then binned to have at least 50 counts in each energy bin. 

The Reflection Grating Spectrometer (RGS; den Herder et al. 2001) data have been reduced 
using the standard SAS task rgsproc and the most recent calibration files; the total exposure 
times are ~ 97 ksec for both RGSl and RGS2. The RGSl and RGS2 spectra were binned 
at the resolution of the instrument (AA ~ 0. lA). 


2 . 2 . Chandra 

Chandra observed MCG-5-23-16 with the ACIS-S with two relatively short exposures for 
a total of 50 ks (see Table 1). For this study the Chandra observations were made with the 
High-Energy Transmission Grating (or HETG - Markert et al. 1994) in the focal plane of the 
High Resolution Mirror Assembly. The Chandra HETG consists of two grating assemblies, 
a High-Energy Grating (HEG) and a Medium- Energy Grating (MEG); the HEG affords the 
best spectral resolution in the ~ 6 — 7 keV Fe-K band currently available (~ 39 eV, or 
1860 km s _1 FWHM at 6.4 keV). The MEG spectral resolution is only half that of the HEG. 
The HEG also has higher effective area in the Fe-K band. The HEG and MEG energy bands 
are ~ 0.9 — 10 keV and ~ 0.4 — 8 keV respectively, though the effective area falls off rapidly 
with energy near both ends of each bandpass. The Chandra data were reprocessed with CIAO 
version 3.2.1 1 and t CALDB version 3.0.1. Spectral redistribution matrices (rmf files) were 
made with the CIAO tool mkgrmf for each arm (—1 and +1) for the first order data of each of 
the gratings, HEG and MEG. Telescope effective area files were made with the CIAO script 
f ullgarf which drives the CIAO tool mkgarf . Again, separate files were made for each arm 
for each grating for the first order. The effective areas were corrected for the time-dependent 
low-energy degradation of the ACIS CCDs using the option available in the mkgarf tool in 
the stated version of the CIAO and CALDB distribution. Events were extracted from the 
— 1 and +1 arms of the HEG and MEG using strips of width ±3.6 arcseconds in the cross- 
dispersion direction. Lightcurves and spectra were made from these events and the spectral 
fitting described later was performed on first-order spectra combined from the —1 and +1 
orders (using response files combined with appropriate weighting), but keeping the HEG 
and MEG spectra separate. Background was not subtracted as it is negligible in the energy 
ranges of interest. Examination of the image of the entire detector and cross-dispersion 
profiles confirmed that there were no nearby sources contaminating the data. 


1 http://cxc. harvard.edu/ciao 
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In the following, unless otherwise stated, fit parameters are quoted in the rest frame of the 
source and errors are at the 90% confidence level for one interesting parameter (Ay 2 = 2.71). 
Abundances were set those of Wilms et al. (2000). 


3. Spectral analysis 

3.1. The XMM -Newton 0.3—10 keV continuum 

To characterize the X-ray continuum of MCG-5-23-16 we first fitted the 2-10 keV pn 
data, with a simple redshifted power-law model, modified by Galactic (Nh = 8 x 10 20 cm -2 ; 
Dickey & Lockman 1990) and local absorption. For this first fit, we fixed the photon index 
(r) to 1.8 and we ignored the 5-7.5 keV band, where the Fe Ka emission is expected. This 
model is a poor description (y 2 /dof= 1312.5/988) of the X-ray emission and its extrapolation 
leaves strong residuals in the soft (E< 1 keV) band and, as expected, at the energy of the Fe 
Ka emission line (see Fig. 1). Furthermore this fit can not account for the pn emission above 
8 keV, which is indicative of a harder X-ray continuum or it could be due to the presence of 
emission due to Compton reflection. Leaving the photon index free to vary we find T ~ 1.65 
and Nh ~ 1-3 x 10 22 cm“ 2 . This is a more acceptable fit (y 2 /dof= 1066.3/987), but it can 
only account for the excess above 8 keV, while it cannot model the Fe Ka line and the soft 
emission. Indeed the extrapolation of this model to the 0.5-10 keV band does not provide 
an acceptable fit (y 2 /dof= 6024.2/1794). 

In order to model the soft X-ray emission, we added to this model a soft power-law compo- 
nent absorbed only by the Galactic column. This model is still too simple to describe the 
overall emission of MCG-5-23-16, however it shows that the soft X-ray emission could be 
explained with scattering of the primary power-law continuum. The photon index of this 
soft component is found to be steep, T = 3.13 ± 0.10, and even at the CCD resolution of the 
pn instrument the power-law model leaves line like residuals (see Fig 2; black data points). 
In particular an emission line is required by the data (Ay 2 = —37) at 0.92 ± 0.02 keV with 
a flux of 2.1 x 10 _5 photons cm -2 s -1 . This emission line and the steep power-law continuum 
is confirmed by the MOS1 and MOS2 data (see Fig. 2, red and green data points) and also 
by the simultaneous Suzaku observation (Reeves et al. 2007). This model gives a 2-10 keV 
observed flux of ~ 8.0 x 10~ n erg cm -2 s"" 1 and an observed luminosity of ~ 1.3 x 10 43 erg s -1 . 
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3.2. The RGS spectra: soft X-ray spectrum dominated by emission lines 

In order to investigate if the apparent steep soft X-ray photon-index could be due to the 
presence of emission lines which are unresolved at the pn CCD resolution, we analyzed the 
RGS data. Indeed, thanks to the long exposure (~ 100 ksec), the RGS1 and RGS2 spectra 
have enough statistics to attempt a spectral analysis (a total of ~ 17500 net counts between 
RGS1 and RGS2). The first inspection of the RGS data reveals the presence of several soft 
X-ray emission lines as well as the energy cut-off at ~ 1 keV due to the rest frame absorption. 
We then rebinned the RGS spectra in constant wavelength bins at the spectral resolution 
of the instrument (AA ~ 0.1 A) and we use the C-statistic (Cash 1979) available in xspec 2 
for the spectral fit, because with this choice of binning, we have some bins with less than 20 
counts. 

We first fitted the RGS spectra with the baseline model obtained with the pn spectrum. 
This model consists of two components: a primary absorbed power law and a scattered soft 
power law component absorbed only by Galactic absorption; both photon indices have been 
fixed to the value found for the primary power law component (r = 1.82, see Section 3.3). 
Overall, this model provides a reasonable description of the RGS continuum, however line- 
like residuals are present below 1 keV. 

We then kept the AGN baseline model parameterized with the above model and we tested 
two different additional components for the emission below 1 keV, which are: (a) a multi- 
temperature thermal emission model with variable abundances for different elements (Mewe 
et al. 1985) or (6) several unresolved emission lines. The first model represents the emission 
expected from a collisionally ionized plasma; the latter reproduces the emission due to ma- 
terial photoionized by the AGN. 

When modeling with the multi temperature model, we found that the data could be fitted 
with two thermal components with kT\ = 0. 441^4 keV and kT^ = 0.151q'o 7 keV (AG = 52 
with respect to the absorbed power law model, for 9 additional parameters). The flux and 
luminosity of these thermal components are F( 0 . 5 _ 2 )kev ~ 9.2 x 10 -14 erg cm -2 s -1 and 
L(o.5-2)keV ~ 1.5 x 10 40 erg s -1 , which are consistent with the possible emission from the 
galaxy. The data below 1 keV still require a contribution from the scattered power law 
component of the AGN baseline model, which has F(o. 5 - 2 )keV ~ 1.4 x 10 -13 erg cm -2 s _1 ; 
allowing the photon index of this soft power law to vary we found that its photon index is 
no longer unusually steep (T ~ 2.3). 

Although not well constrained, the abundances required with this model are found to be 
low, Z ~ 0.2 Zq. In particular iron is found to be underabundant, with only an upper limit 
of Q.2Z q . This is due to the lack of Fe-L shell emission lines and is at odds with the flux 


2 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/manuai/manual.html 
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measured for Fe Ka emission line (see section 3.3 and Table 3). Furthermore this value is 
also in contrast with the Fe abundance measured with the neutral iron edge from the reflec- 
tion component in the simultaneous Suzaku data ( Zp e = 0.4 ± 0.1 Zq at the 90% confidence 
level; Reeves et al. 2007). 

Instead we fitted the RGS spectra adding to the AGN baseline continuum model several 
unresolved emission lines, fixing both the soft and the hard power-law photon indices to 
the value found for the AGN primary power- law component (T = 1.82 see section 3.3). To 
account for the excess of counts below 1 keV, five lines are required (AC = 70); their fluxes 
are listed in Table 2, together with their EWs which range from ~13 eV (N vil) to ~ 57 eV 
(Ovil). The most likely identifications are 2— >1 emission lines from H and He-like O, Ne 
and N (see Fig. 3). 

We tentatively allowed the lines’ width to vary and found the O vil Hea, the O vil RRC, 
and the Ne IX Hea lines to be marginally resolved, while the N VII Lya and the O vm Lya 
lines are unresolved. Although a quantitative measure is beyond the statistics of the present 
data, from the width of the Ovil RRC feature (E> 0.739 keV) we derive an upper limit 
of kT < 24 eV on the recombining electron temperature. It is worth noting that this low 
value indicates that the soft X-ray emission originates in a photoionized rather than in a 
collisionally ionized plasma (Liedahl & Paerels 1996; Liedahl 1999). 

The O VII Hea and the Ne IX Hect lines are both triplets, which with the present statistics, 
can not be resolved into their forbidden and resonant components. However it is worth 
noting that for both the O vil Hea and the Neix Hea lines the energy centroids are close to 
the energy of the rest- frame forbidden line (see Table 2). This suggests a major contribution 
from forbidden lines in both of the triplets with respect to the resonant lines. 

In the case of Ne IX Hea line we cannot derive any quantitative information on the width 
of this triplet, since it is close to the Fexix 3d-2p blend (0.913-0.926 keV) and to further 
complicate the analysis its energy is close to the low energy photoelectric cut-off present in 
the MCG-5-23-16 spectrum due to the local absorption (Nh ~ 10 22 cm~ 2 ). 

In the case of the O VII Hea we measured a width of a = 4.6lf ieV, which is probably due to 
the presence of the weak resonant component. In order to confirm that the detected emission 
is dominated by the forbidden line, we then added a second line and kept the line energies 
frozen for the forbidden (0.561 keV) and the resonant (0.574 keV) lines. The measured ratio 
between the flux of the forbidden and resonant lines (for the latter we use the 90% upper 
limit) is £ 2.55; which again is an evidence of a strong contribution from a photoionized 
plasma (Porquet & Dubau 2000). Furthermore with this model we found that the width of 
the forbidden line is now unresolved (a < 4.6 eV). 

Finally an inspection of the values of the centroid energies of the detected lines (see Table 2) 
shows that there are no strong shifts between the theoretical and observed values. The mea- 
sured shifts of these lines give a value of 1 eV (which corresponds to a velocity 2000 
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km/s). 

In order to test whether these lines can be explained with emission from optically thin gas 
photoionized by the AGN we replaced the unresolved emission lines with a grid of photoion- 
ized emission models generated by xstar (Bautista &; Kallman 2001), which assumes a T ~ 2 
illuminating continuum and a turbulence velocity of cr v = 100 km/s. We found that the RGS 
data are well explained with this model with an ionization parameter log £ = 1.271q 13- We 
then let the photon index of the scattered power- law vary and we found that, although it 
is not well constrained, the value is now similar (r = 1.89i°;g|) to the primary AGN power 
law emission, in agreement with the scattering hypothesis. Furthermore, the abundance of 
photoionized iron (~ 0.4 Z & ) is now similar to the value obtained with the detailed fitting 
of the Fe K a line and the Compton reflection hump detected with the Suzaku observation 
(Reeves et al. 2007). 

As a final test, we applied this best fit model to the pn, MOS1 and MOS2 data. We kept the 
abundances fixed to the values measured with the RGS spectra and we let only the normal- 
ization and the photon index vary. This model is now a good description of the soft spectrum 
and no strong residuals are present. Finally with this model we found that the scattered 
component has a T = 2. 24^24 an d F( 0 .5-2)keV ~ 4 x 10~ 13 erg cm -2 s -1 . This corresponds 
to 0.5% of the un-absorbed flux of the primary AGN component and the emission due to 
the photoionized gas (L(o.5_2)k e v ~ 10 40 erg s -1 ) is 0.1% of the AGN emission. 

To conclude, although from a statistical point of view either the multi-temperature thermal 
emission model or the photoionized plasma model (xstar) give similar results, different 
diagnostics suggest that, as already seen in other Seyfert2 galaxies (Bianchi et al. 2006; 
Guainazzi & Bianchi 2007; Iwasawa et al. 2003) the soft X-ray emission is likely due to pho- 
toionized plasma which could be associated with the Narrow Line Region. The properties 
of the soft X-ray emission of MCG-5-23-16 that favor this interpretation are: the lack of 
a strong Fe-L shell emission, the detection of a narrow RRC feature from O VIII and the 
stronger contribution from the forbidden line in the O VII Hea triplet. 

In order to assess the extent of the soft flux, we used the Chandra observations and created 
an image combining the 0. 3-1.0 keV photons from the zero-order of both observations. We 
created a point-source function (PSF) using the MARX 3 Chandra simulator. We then fit a 
model to the Chandra image consisting of a constant component, to account for the back- 
ground, and two Gaussian components with the centroid positions tied together, to account 
for both nuclear and extranuclear emission. The model was convolved with the PSF and 
then compared to the data using the Cash statistic. Initially the image was fit with cr x = a y , 
i.e., circular Gaussian models. This gave o = 0.39'%o'i5 for the nuclear component and 


3 http://space. mit.edu/CXC/MARX/ 
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a = 1.5 7" to! 3 i f° r the extranuclear component. The count rates for the two components 
where 1.4 x 10 -3 counts s -1 and 1.2 x 10 -3 counts s -1 . Note that some fraction of these 
extents is likely due to residual error in the aspect solution which would effectively increase 
the PSF. Unfortunately there are no other on-axis point sources bright enough in the field 
to check this. Allowing the extra-nuclear component to be elliptical did not improve the 
fit significantly, however it reduced C by 7.2 for 2 additional parameters (the additional a 
parameter and the rotation angle), which is significant at the 2 a level. It is likely, however, 
that this asymmetry may be due to aspect errors. As a consistency check we fitted the 
zero-order image from the 2000 Chandra observation with the same model. In this case al- 
lowing for ellipticity in either component did not improve the fit significantly. In this fit the 
best-fitting parameters were a = 0.56"io!i2 f° r the nuclear component and o = 2.2"^^ for 
the extra-nuclear component, which is consistent within the errors. This is indicative that 
around half of the soft X-ray emission is due to the central point like source and half to an 
extended component. Assuming the current cosmology (Ho = 71 km s -1 Mpc -1 , f7\=0.73, 
and G m =0.27) the soft X-ray emission of this latter component originates within ~ 0.7 kpc; 
this value is in agreement with a possible association with the NLR; indeed this extension 
is comparable to the extension of the [O ill] A5007 derived with HST data (Ferruit et al. 2000). 


3.3. The iron K band 

We then considered the hard X-ray emission of MCG-5-23-16, using the dual power law 
continuum as described above and only the Gaussian emission line at ~ 0.9 keV. 

In Fig. 4 we show the residuals left by the absorbed power law model (with T = 1.65) at the 
energy of Fe K band. These residuals clearly reveal the presence of a strong narrow core at 
the expected energy of the Fe Ka (6.4 keV) and broad wings, which extend from ~ 5.7 keV 
to ~ 7 keV. The pn data also show a narrow emission line at ~ 7 keV, due to Fe K/3 and 
a drop at 7.1 keV probably due to a reflection edge. The presence of this latter component 
was already suggested with the previous short X.MM- Newton observations (Dewangan et 
al. 2003), however the short exposure of these observations together with the lack of any 
simultaneous observation above 10 keV did not allow the authors to put strong constraints 
on this feature. 

We adopted the best-fit model obtained by Reeves et al. (2007) from the simultaneous Suzaku 
observation for the underlying continuum in order to derive the Fe Ka emission line prop- 
erties. Indeed Suzaku’s broad band energy range (0.4-100 keV) allowed these authors to 
measure the amount of Compton reflection and thus better constrain the continuum in the 
2-10 keV energy band. This model is composed of: a primary absorbed power-law com- 
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ponent with an exponential cut-off at high energies (> 200 keV) and a component due to 
reflection from neutral material (the PEXRAV model in Xspec, Magdziarz & Zdziarski 1995). 
The parameters of this reflection component are: a reflection fraction, which is defined by 
the subtending solid angle of the reflector R = Q/ 2n = 1.1; an inclination angle i = 45° and 
abundance Z = OAZq. 

When fitting this model to the pn data we keep the values of Z, R and the cut-off energy 
fixed, since they cannot be determined using the lower energy band pass of XMM-Newton. 
After including the reflection component the residuals no longer show the deep edge at 7.1 
keV which is well modeled with the reflection component (see Fig. 5; upper panel). With this 
model we found that the primary power law component has a photon index T = 1.82 ± 0.01, 
absorbed by a neutral column density of Nh = 1.49 ± 0.01 x 10 22 cm -2 . 

To model the Fe line we first added narrow Gaussian lines at the energies of Fe Ka and Fe 
K (3. For this latter line we kept the energy fixed at 7.06 keV and tied its flux to be 12% 
of the Fe Ka flux. This model clearly leaves an excess of counts at the energy of the 6.4 
keV Fe Ka line (see Fig. 5; middle panel), which can be accounted for by including a broad 
Gaussian line or a relativistic disk-line component (see section 3.3.2). 


3.3.1. Chandra Observation of the Narrow Core. 

In order to measure the parameters (i.e. strength and profile) of the broad compo- 
nent we first derived the width and flux of the narrow core using the simultaneous Chandra 
observations. We combined the Chandra ±1 MEG and HEG 1st order spectra of both ob- 
servations. The combined spectra were rebinned at the maximum spectral resolution of the 
instruments (AA = 0.012 A and 0.023 A for HEG and MEG respectively) and the spectral 
fit was minimized with the C-statistic (Cash 1979). 

We then adopted the Suzaku best fit model for the underlying 2-8 keV continuum. Thanks 
to the high resolution of the MEG and HEG instruments, the Chandra HETG residuals 
clearly reveal the presence of a narrow core at 6.4 keV (see Fig. 6), which is best modeled 
with a Gaussian line at 6. 40^01 keV and EW~ 80 eV (AC = 43). With this model we 
measured a width of a — 32ljg eV, which corresponds to a velocity width of cr v ~ 1400 km/s 
and is thus indicative of a possible origin from the molecular torus. 

Taking into account that the measured width of the Fe Ka narrow core could be due to the 
presence of the broad component, we added a second Gaussian line. Although the fit does 
not statistically improve, we found an EW of ~ 60 eV (flux ~ 5.6 x 10 -5 photons cm -2 
s -1 ; a ~ 0.4) and ~ 70 eV (flux ~ 5.8 x 10~ 5 photons cm -2 s _1 ) for the broad and nar- 
row components respectively, which are in good agreement with the values measured with 
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the XMM-Newton (see Table 3) and Suzaku spectra (Reeves et al. 2007). With this model 
the narrow core is no longer resolved and we can place only an upper limit on the width 
of a < 43 eV, which correspond to a FWHM < 5000 km/s (at the 90% confidence level) 
in good agreement with the upper limit measured with a previous Chandra observation of 
MCG-5-23-16 (Balestra et al. 2004). 


3.3.2. The broad Fe line 

We then adopted the Chandra upper limit on the width of the narrow core for the 
XMM-Newton fits. The broad component was first modeled adding a second Gaussian line. 
We also added a Compton shoulder at 6.3 keV, with its normalization set to 20% of the Fe 
Ka flux (Matt 2002). The fit improves (A% 2 = 44 for 3 additional parameters; x 2 = 1932 
for 1785 dof). The broad Gaussian component (E=6.22 /q;]; 6 keV) has an EW of 66/}g e V 
and a width of a = 0.42/ojo keV; which corresponds to a velocity of a v ~ 20000 km/s 
(FWHM~ 40000 km/s). 

We then tested a relativistic diskline model ( DISKLINE in XSPEC; Fabian et al. 1989); 
this code models a line profile from an accretion disk around a Schwarzschild black hole. 
The main parameters of this model axe the inner and outer radii of the emitting region on 
the disk, and its inclination. The disk radial emissivity is assumed to be a power-law, in 
the form of r~ q . For the fit we fixed the outer radius to be 400R g (with R g = GM/c 2 ) and 
the emissivity to be q — 3. Finally we assumed the line to be from neutral Fe Ka. From a 
statistical point of view this fit gives a similar result to the model with a broad Gaussian line 
(x 2 = 1928 for 1785 dof), however the high velocity inferred from the width of the Gaussian 
line is indicative that the line must be produced close to the central black hole; i.e. within 
100Rg and thus inside the Broad Line Region. 

With this model we found that the inner radius is Ri n = 48/20 Rg, an< ^ the inclination angle 
is i = 41°/^; while the EWs of the broad and narrow components are EWDiskiine = 53/^ 
eV and EWn = 64 ± 6 eV. If the constraint on the disk emissivity is relaxed and a flatter 
emissivity is assumed ( q = 2) then a disk inner radius of 6R g is allowed by the present data. 
The ratio between the data and this best fit model is shown in Fig. 5 (lower panel), an ab- 
sorption line near 8 keV is the only residual. Upon adding an absorption line to our best fit 
model improves the fit (A% 2 = 20 for 2 additional parameters, \EW\ ~ 30 eV, E~ 7.9 keV). 
The parameters derived for the diskline do not change significantly; the main difference is 
a slight reduction in its EW, which is now 46/ eV. Taking into account that this absorp- 
tion feature is indicative of the presence of an ionized absorber (see section 4.1 and 5.3), 
we tested if the presence of a more complex absorber can mimic the profile of the detected 
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broad component. We found that the inclusion of a two layers of absorption, characterized 
by a high (log£ = 3.7) and lower (log£ = 2) ionization level, does not impact the detection 
of the broad component and its parameters. 


4. Variability of the Iron line and continuum 

During the XMM-Newton observation, the 2-10 keV flux of MCG-5-23-16 varied from 
~ 7 x 10~ n erg cm -2 s -1 to ~ 9 x 10~ n erg cm -2 s -1 . In order to investigate the possible 
variability of the line properties and continuum shape, we first tested if there was any clear 
difference between the spectra extracted when the source was in a higher and lower flux state. 
We extracted spectra using a 2-10 keV threshold of < 4.6 cts/s (F( 2 _iokeV) ~ 7.4 x 10 -11 
erg cm -2 s -1 ) and > 5.5 cts/s (F( 2 _iokeV) ~ 8.8 x 10~ n erg cm -2 s -1 ). We then fitted both 
spectra with the previous best-fit model but replacing the diskline with a broad Gaussian. 
We found no evidence of variability of the broad or of the narrow component, indeed their 
normalizations are consistent within the errors (see Table 3). Furthermore, the width of the 
broad line is constant. In order to confirm that the line is not strongly variable we checked 
the difference spectrum, obtained by subtracting the low from the high state data. The 
difference spectrum can be essentially modeled with an absorbed power law having a photon 
index T = 1.80 ± 0.09, no strong residuals are left in the soft band as well as in the iron K 
band. 

As a second test to assess the possible variability of the Fe emission complex, we divided the 
XMM-Newton observation into 5 intervals with a duration of 20 ksec each. 

We tied all model parameters except the normalization of the primary power law component: 
the model gives a statistically acceptable fit for all the 5 spectra. Fig. 7 shows the 5. 5-8. 5 
keV residuals to this model for all five intervals; no strong deviations from the model are 
required at the energy of the iron Kct line. We then allowed the normalization of the narrow 
and broad components free to vary. Fig. 8 (panel a and panel b) shows the flux variability of 
both the narrow and broad component. There is no evidence of variability during the present 
observation, furthermore the fluxes of both components are consistent within the errors with 
the values measured in the previous XMM-Newton and Chandra observations (Dewangan et 
al. 2003; Balestra et al. 2004). This lack of variability can be explained taking into account 
that MCG-5-23-16 is not highly variable on both short and relatively long time-scales, indeed 
the source remained at a similar flux level (7-9 xl0 _11 erg cm -2 s -1 ) for the last 10 years. 
The most striking feature that appears to be variable during the XMM-Newton observation 
is a possible absorption line at ~ 7.7 keV (E~ 7.66 keV, observer frame). This feature 
is prese nt in Ihe average EPIC pn spectrum but as shown in Fig. 7, it is strongest in the 
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third spectrum (40 ksec after the beginning of the XMM- Newton observation). To illustrate 
this in Fig. 8 (panel c) we compare the intensity of this absorption feature during the five 
intervals. For this purpose we modeled the absorption with an inverted Gaussian line and 
fixed the energy to the best-fit value found with the spectral analysis of the third interval 
(E~7.7 keV, |/| ~ 3.2 x 10 _5 photons cm _2 s _1 , EW= 52 ± 15 eV). The line is clearly variable 
and it appears to be stronger during the interval with the slightly higher 2-10 keV flux (see 
Fig. 8 panel d) while it is barely detected in the other four spectra. 


4.1. A variable absorption feature at 7.7 keV 

4-1.1. Epic-pn background and calibration checks. 

Before attempting any further modeling, we performed several tests to exclude that the 
7.7 keV absorption feature is due to inappropriate background, binning or pattern selection. 
The EPIC-pn background near this energy range presents two instrumental lines due to Cu 
(8.05 keV) and Ni (7.48 keV) Ko emission lines and an inadequate background selection 
could in principle cause spurious absorption features. However several arguments exclude 
this possibility. First of all the net count rate (~ 2.56 cts/s) of MCG-5-23-16 in the 5-10 keV 
band is ~ 300 times greater than the background (~ 8 x 10 -3 cts/s). Second, the sporadic 
nature of the feature is indicative that the feature can not be an artifact of the background or 
calibration of the EPIC-pn. Finally there was no background flaring activity during this time 
interval. We conclude that the feature is not due to instrumental or external background. 

In order to exclude the possibility that the 7.7 keV feature is due to a binning effect we 
rebinned the pn data of the third interval with constant energy binning of 80 eV. This choice 
corresponds to about half of the energy resolution EPIC-pn camera in this energy range 
(FWHM~ 150 eV at 6.4 keV; see the XMM-Newton Users’ Handbook, Ehle et al. 2006 4 ). 
As shown in Fig. 9, the residuals left by the time averaged best fit model confirm the presence 
of the feature, thus excluding the possibility that it is an artifact of the choice of the binning. 
To exclude a pattern selection effect we then compared the pn spectra extracted with the 
pattern 0-4 and pattern 0 selection criteria. Though the latter has 30% fewer counts we 
found no significant difference in the absorption line parameters (Ay; 2 = 26 for 2 dof; E~ 7.7 
keV, EW~ 60 eV). Finally the presence of an absorption feature is confirmed by the MOSl 
and MOS2 spectra extracted in the same time interval. Although, due to the lower photon 
statistic, the significance of the absorption line is lower in these spectra, both the flux and the 
energy of the feature are consistent (|/| = 2.2 ± 1.8 x 10~ 5 photons cm _2 s _1 E=7.4±0.2 keV 


4 Vit.t.pY/vmm vilspa. p.sa. ps/ftvtfirna.l/ymm user snppnrt./Hnr.nmpnt.a.t.inn/nhh/YMM TTTTR ht.ml 
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Ay 2 = 6; see Table .4) with the values found with the pn data 5 . Furthermore the presence of 
the absorption feature is confirmed by the simultaneous Suzaku observation, indeed a weak 
absorption feature is present in the time averaged spectrum (see panel c of Fig. 6 in Reeves 
et al. 2007). The absorption line is weaker than in the XMM observation, which could be 
explained if we take into account the apparent sporadic nature of the feature, with a dilution 
effect due to the longer duration time (~220 ksec) of the Suzaku observation in the Earth 
orbit. Although the absorption line is not well constrained in the time-averaged Suzaku 
spectrum, the energy of the line at 7.8 ± 0.1 keV is coincident with the XMM-Newton data, 
while the flux of the line is weaker at —1.8 ± 0.9 x 10 -5 photons cm~ 2 s -1 (see Table 4). 


4-1.2. Modeling the absorption feature of the third segment. 

We first fitted the absorption feature adding a Gaussian shaped absorption line keeping 
its width fixed to a = 0.1 keV. The addition of this line improves the fit with a Ay 2 = 33 
for 2 additional parameters (y 2 /do/= 1451/1505). The line energy is 7.72 ± 0.06 keV with a 
|EW| of 52 ± 15 eV (|/| = 3.2 ± 0.9 x 10~ 5 photons cm -2 s -1 ). Leaving the width of the line 
free does not improve the fit significantly (Ay 2 = 5 for one additional parameter). With this 
fit we find <7 = 0.2 ±0.1 keV, EW=78 ± 29 eV and an energy consistent with the previous 
best fit (E=7.71 ± 0.08 keV). In Fig. 10 we show the confidence contour plot of the line 
parameters (rest frame energy and intensity) with the line width left free to vary. 

We also attempted to fit the absorption feature replacing the Gaussian line with an edge 
due to K-shell absorption from partially ionized iron. This model gives a best-fit energy of 
7-33lo 22 and an optical depth of r = 0.09 ±0.03. The fit is statistically acceptable; however, 
it is worse than the Gaussian absorption model (y 2 /do/=1465/1505, which correspond to a 
Ay 2 = 14 worse compared to the absorption line) and more importantly, it is unsuccessful 
at modeling the residuals at 7.7 keV. 


5 Unfortunately the Chandra observations do not overlap with this segment of the XMM-Newton spectrum. 
However two possible weak absorption feature appear to be present at the rest frame energy of about 7.3 keV 
and 7.4 keV (see Fig. 6) suggesting possible variability of the absorber, although the statistical significance 
of these features is low. 
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4-1-3. The significance of the detection of the absorption line 

By applying the standard two-parameter F-test to the drop in y 2 of 33 for the addition 
of an inverted Gaussian at 7.71 keV, we found a null hypothesis probability for adding this 
extra component of ~ 4 x 10 -8 . However, as discussed by Protassov et al. (2002) the F-test 
applied in this way could overestimate the true significance of the detected absorption line. 
In particular the F-test does not take into account the number of time bins in which the line 
is searched as well as the range of energy where the line might be expected (see Porquet et al. 
2004). To asses the significance of the detection we then performed Monte Carlo simulations 
as described in Porquet et al. (2004) and in Markowitz et al. (2006), for a similar case of a 
blue-shifted Fe Kcc absorption line detected in the XMM -Newton observation of IC 4329a. 
We assumed as our null hypothesis model the best fit model with no absorption feature, and 
we simulated 1000 spectra with the photon statistics expected for a 20 ksec exposure. Each 
simulated spectrum was then fitted with the null hypothesis model to obtain a y 2 value and 
systematically searched for an absorption line over the 4-9 keV energy range with a step of 
0.1 keV in the centroid energy of the possible line and re-fitted at each 0.1 keV step. We 
then obtained for each simulated spectrum a minimum y 2 and created a distribution of 1000 
simulated values of the Ay 2 (compared to the null hypothesis model), which were used to 
construct a cumulative frequency distribution of the Ay 2 expected for a blind line search 
in the 4 — 9 keV range. Not a single fake spectrum had a | Ay 2 1 > 33, thus the inferred 
probability that the null hypothesis model was correct is < 0.1%. Taking into account the 
number of intervals (5) in which the observation has been split into we derive that the line 
detection is significant at > 99.5%. 

Finally, we performed a similar Monte Carlo simulation to test the statistical significance of 
the absorption feature in the mos spectra. For simplicity we ran the simulation on MOS1 
only and we found that in this case the significance is only ~ 61%, which mainly due to the 
lower S/N of the mos spectra at this energy. However the fact that the absorption feature 
is detected by the pn, both mos and Suzaku suggests that the feature is likely real and not 
an artefact. 


4-1-4- The ionized absorber model. 

As already discussed from X-ray spectroscopic observation of several other Seyfert galax- 
ies (i.e. MCG-06-30-15, Young et al. 2005; NGC 3783, Reeves et al. 2004; Mrk 509, Dadina et 
al. 2005; Mrk 335 Longinotti et al. 2007; E1821+643 Yaqoob & Serlemitsos 2005; IC 4329a 
Markowitz et al. 2006) and QSOs (i.e. PG 1211+143, Pounds et al. 2003; PDS456 Reeves et 
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al. 2003; APM 08279+5255 Chartas et al. 2002; PG1115+080 Chartas et al. 2003) a likely 
candidate for the 7.7 keV feature is blue- shifted K-shell absorption due to He- or H-like iron. 
In particular if we assume that the line is due to H-like iron (Fexxvi Lya at 6.97 keV) the 
observed blue-shift suggests that the absorber is outflowing with a velocity of the order of 
0.1c. 

In order to obtain a more physical representation of the absorber we replaced the Gaus- 
sian line with a model comprising a grid of photoionized absorbers generated by the XSTAR 
photoionization code (Bautista & Kallman 2001). For the absorber we assumed a one zone 
photoionization model with half solar abundances and a turbulence velocity of 1000 km s _1 . 
The free parameters of this model are: the column density (Nh), the outflowing velocity of 
the absorber (u out ) and the ionization parameter (£ = L/nr 2 ). To reproduce the absorp- 
tion feature a column density of ~ 8 x 10 22 cm “ 2 and an ionization state characterized by 
log£ = 3 . 7 I 0 J erg cm s -1 are required with an outflow velocity of 0.09 ± 0.01c (~ 30000 
km/s). A plot of this best-fit model is shown in Fig. 11 (model A). The column density 
is not well constrained and we can place only a lower limit of > 2 x 10 22 cm ~ 2 . At this 
ionization level the Fe K-shell absorption is indeed mainly due to Fe XXVI and is consistent 
with absorption from highly ionized iron outflowing at ~ 0.1 c with respect to systemic. 

A lower velocity outflow could in principle be obtained assuming that the feature is due 
to a Fe xxiv ls-3p line at 7.78keV. This is illustrated in figure 11 (model B), which illus- 
trates an XSTAR model with a column density of 7V H = 10 23 cm~ 2 , an ionization parameter 
of log£ = 3.0 erg cm s _1 and no velocity shift. The absorption line at 7.8 keV indeed cor- 
responds to the above Fe xxiv 1 — > 3 transition. However at this lower ionization we would 
expect to also see the stronger absorption due to the Fe xxiv ls-2p absorption line at 6.67 
keV, which is not detected. Indeed as illustrated in the figure, we would expect to detect a 
strong absorption trough due to a blend of the 1 — > 2 transitions of Fe xvn-xxv, between 
6.5-6.7 keV, which is not observed in the MCG -5-23-16 spectrum. Furthermore in this ion- 
ization regime, several strong absorption lines from iron L-shell (2 — » 3) transitions as well 
as He and H-like Si/S K-shell lines are expected near 1-2 keV, which are not observed in 
either the XMM-Newton, Suzaku or Chandra HETG spectra. Therefore this lower velocity 
solution appears to be ruled out. 

Finally it is possible to have a low ionization K/3 absorption line without strong Ka from 
the same species, when the Fe L-shell is fully occupied. This scenario is shown in Figure 11 
(model C) for an absorber with a column density of N h = 10 23 cm -2 , an ionization parame- 
ter of logl; = 1.5 erg cm s -1 and no velocity shift. There is no Is - 2 p absorption between 
6.4-7.0keV as the L-shell is filled, however there is an absorption line from 7.1-7.2keV due 
to 1 — > 3 transitions from Fe< xvn. In this scenario a blueshift of ~ 0.08c would still be 
required to model the absorption line at 7.7keV in the spectrum. Furthermore the column 
density actually required to mod e l the EW of the Kfi absorption feature is Nh >> 10 2 3 cm~ 2 , 
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which would introduce too much continuum (bound-free) absorption below 6 keV, inconsis- 
tent with the pn data. 

Therefore we conclude that the fast (0.1c) high ionization outflow is the most likely model 
to account for the absorption feature at 7.7 keV. Moreover, when we compare the 2-10 keV 
continuum and the neutral Nh measured during the 3rd interval with the other intervals and 
with the average spectrum, we do not find any statistically significant difference (ANjj < 10 21 
cm -2 ), which rules out the possible presence of variable neutral and/or low ionization ab- 
sorber. 


5. Discussion and Conclusion. 

We have presented the results from XMM- Newton and Chandra observations of MCG- 
5-23-16, which are part of a simultaneous campaign conducted in December 2005 also com- 
prising Suzaku and RXTE. The 0.5-10 keV continuum of MCG-5-23-16 resembles at first 
order the canonical X-ray emission expected from a Compton thin Seyfert 2 galaxy: an ab- 
sorbed (Nh = 1-5 x 10 22 cm -2 ) power law component (T = 1.82), which emerges at energy 
XL 1 keV, and a steep soft excess, which is well fitted by a power law component plus several 
emission lines from O, Ne and N. The XMM -Newton observation of MCG-5-23-16 confirms 
the presence of the Fe Ka emission complex, which is well described by a narrow Fe Ka 
emission line superimposed on a relativistically broadened component. The simultaneous 
Suzaku observation provided us an accurate description of the underlying continuum which 
allowed us to perform detailed modeling of the Fe K emission line complex. Finally due to 
the sufficiently long duration, the XMM-Newton-RGS spectra have enough photon statistics 
to investigate the origin of the soft X-ray emission. 


5.1. The origin of the soft X-ray emission 

The analysis of the EPIC-pn and MOS spectra of MCG-5-23-16 revealed the presence, 
below 1 keV, of a soft excess with respect to the primary nuclear emission. This soft ex- 
cess can be well fitted adding an un-absorbed power law component to the primary AGN 
emission, the photon index of this power law is found to be steeper (r ~ 3.1) than the 
primary AGN component (r ~ 1.8), and even at the EPIC CCD resolution an emission line 
is detected around 0.9 keV. A soft excess below 1 keV is not unusual in obscured Seyfert 
galaxies like MCG-5-23-16 (see Bianchi et al. 2006 and references therein) and it has been 
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already suggested that it could be due to a superimposition of scattered emission into the 
line of sight by ionized gas plus several emission lines from highly ionized (He and H like) 
elements (i.e Mrk 3: Sako et al. 2000; Pounds & Page 2005; Bianchi et al. 2005; Circinus: 
Sambruna et al. 2001, NGC 1068: Kinkhabwala et al. 2002; Ogle et al. 2003; Brinkman et al. 
2002; NGC 4507: Matt et al. 2004). Key diagnostics to understand the origin of this X-ray 
emission when a high resolution spectrum is available are: the detection of RRC transitions, 
the detection of enhanced K-shell emission lines from H-like and He-like ions, the detection 
of Fe L-shell emission and the ratio between the forbidden and the recombination transition 
in the He-like triplets. 

This long XMM -Newton observation provided for the first time RGS spectra for MCG-5-23- 
16 with sufficient photon statistic to perform a detailed modeling of the soft X-ray emission 
and allowed for the first time in this object the detection of the O VII RRC. The width of this 
emission line indicates that the recombining electron temperature is a few eV (kT < 24eV). 
This is suggestive that the emitting plasma is photoionized rather then collisionally ionized 
(Kinkhabwala et al. 2002) and thus it is in agreement with a scenario where the soft X-ray 
emission is due to the AGN activity. Our analysis of the RGS spectra of MCG-5-23-16 con- 
firms the detection of the O VII and O VIII lines previously reported by Guainazzi & Bianchi 
(2007) and the line fluxes axe in agreement with the measurement obtained with this analysis 
of the previous short XMM -Newton observations of MCG-5-23-16. The photon statistics ob- 
tained with the current XMM- Newton observation allowed us also to detect the Neix triplet. 
We cannot exclude on a statistical ground the possible presence of thermal emission; indeed 
the spectra can be equally modeled replacing these emission lines with a multi temperature 
thermal model; which represents the emission due to a collisionally ionized plasma. However 
the “AGN” model (scattered power law component plus several photoionized emission lines) 
is preferred because of the low ion abundance obtained ( Z < 0.2 Z Q ) in the thermal model, 
which reflects the lack of a strong Fe L-shell emission with respect to the fluxes of the Oxygen 
lines. 

We therefore conclude that the most likely the origin of the soft X-ray emission is due to a 
plasma photoionized by the AGN. This plasma must be located outside the Compton-thin 
absorber and, as already suggested for other Seyfert 2 galaxies, it could be coincident with 
the NLR. 


5.2. The Fe Ka emission complex 

This deep XMM -Newton observation of MCG-5-23-16 confirms the presence of broad 
and narrow iron Ka emission linc3, which were reported sinc e the first ASCA observa- 
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tion (Weaver et al. 1997). The Chandra / HETG spectrum clearly reveals a narrow line at 
En = 6.40 ± 0.02 keV with a FWHM< 5000 km/s and a flux of 5.6 ± 0.7 x 10~ 5 photons 
cm -2 s -1 . The intensity of this component is found to be constant, within the errors, during 
this observation and also when comparing with previous observation (< In >= 4.5 x 10 -5 
photons cm -2 s _1 , Balestra et al. 2004; In = 6.5 ± 2.7 x 10 -5 photons cm -2 s -1 , Weaver et 
al. 1997). The constancy in flux of this line together with the limits on the width obtained 
with Chandra are suggestive of an origin from distant matter such as the putative torus. 
Indeed the upper limit on the FWHM corresponds to a distance from the central black hole 
greater then 10 4 R fl . 

The presence of two Compton-thick X-ray reprocessors responsible for the two components 
of the iron line, suggested since the ASCA observation (Weaver et al. 1997) is confirmed 
with this deep XMM- Newton observation and with the deep Suzaku observation (Reeves et 
al. 2007). The geometry inferred for MCG-5-23-16 with this latter observation is discussed 
in detail in Reeves et al. (2007). To summarize, one plausible scenario is that we are seeing 
MCG-5-23-16 through the Compton thin edge of the putative torus, which is Compton thick 
at the plane of the accretion disk. This is in agreement with the inclination ( i = 41°1 2 2) 6 
of the accretion disk derived while modeling the broad line component with the relativistic 
line profile. The column density of 1.5 x 10 22 cm -2 measured with the low energy cut-off is 
thus associated with the thinner absorbing material, e.g. encountered viewing through the 
edge of the torus. We found no evidence of variability of the column density of this ab- 
sorber within this long observation, and also no strong variability is found when comparing 
the column densities measured with the previous observations performed with BeppoSAX 
(Risaliti 2002a), ASCA (Weaver et al. 1997), XMM- Newton and Chandra (Dewangan et al. 
2003; Balestra et al. 2004). This result implies that the absorber is probably far from the 
central black hole and there is no evidence that this absorber is clumpy as suggested for 
other Seyfert 2 (Risaliti et al. 2002b). 

This deep XMM -Newton observation confirms the presence of a relativistically broadened 
iron K<x line; the width derived modeling this component with a Gaussian profile, corre- 
sponds to a FWHM~ 40000 km/s and is suggestive of an origin from the accretion disk. 
The profile of this component is nearly symmetric and can be equally modeled with a broad 
Gaussian or a relativistic profile; in the latter case the derived inner radius is about 20-40i? s . 
Since the advent of XMM- Newton and Chandra one of the most debated issues in the study 
of the broad iron Ka lines has been the fraction of AGNs which clearly show the presence of 
a broad component (Nandra et al. 2006; Guainazzi et al. 2006). Several authors have indeed 


6 Modeling the combined Suzaku and XMM- Newton spectra, Reeves et al. (2007) found an inclination of 
i=53° which is in agreement within the errors with the value obtained with the analysis of the XMM- 
N e.wt.on s pectrum 
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discussed the robustness of the detection in some objects of a broad iron line (i.e. NGC3516 
Turner et al. 2005; NGC3783 Reeves et al. 2004). This controversy emerged when observa- 
tions characterized by high photon statistics showed the ambiguity of modeling the iron Ka 
line when complex absorption is present (Pounds et al. 2003, 2004). Indeed a high column 
density warm absorber can produce curvature in the spectrum at the energy of the iron Ka 
line that mimics the profile of a relativistically broadened emission line. The detection of 
the absorption feature at ~ 7.8 keV shows that a high column density variable absorber (a 
high velocity higly ionized outflow) is present also in MCG-5-23-16 which could in principle 
give rise to ambiguity in the interpretation of the broad component. However, in the case 
of MCG-5-23-16 the availability of a simultaneous observation with Suzaku above 10 keV 
allowed us to tightly constrain the underlying X-ray continuum and to rule out the interpre- 
tation of the broad line as due to unmodeled complex absorption. Indeed the residuals left 
at the energy of the iron Ka, when we take into account the amount of reflection detected 
with Suzaku, cannot be explained by the effect of complex absorption. Furthermore, once 
the absorber responsible for the feature detected at ~ 7.8 keV is accounted for, either in the 
time averaged spectrum or in the third segment of the XMM -Newton observation, a broad 
line is still required by the data with a similar EW and FWHM. Note that the ionization 
parameter of the absorber is required to be high and does not introduce additional spectral 
curvature below 6 keV, which hence does not impact the broad iron line modeling. 

The remaining open questions on the origin of the broad line in MCG-5-23-16 are the rela- 
tively large inner radius derived for the accretion disk, and its lack of variability. The former 
can be explained with several scenarios: the disk could be truncated or missing below 20R S , 
or the inner part of the accretion disk could be so highly ionized that the iron is fully ionized. 
However it worth noting that, as shown by Reeves et al. (2007), assuming a flat emissivity 
(q — 2) a inner radius of about 6R S (in the case of a Schwarzschild black hole) cannot be 
statistically ruled out. 

The latter open issue is the lack of variability of the iron emission line both on short and 
long term time-scales. Indeed the flux of the broad component is found to be consistent 
with being constant both with the short-term time resolved spectroscopy performed within 
this deep observation, and comparing our result with the long term flux history presented in 
Balestra et al. (2004). The strength of the broad component appears to be lower during this 
observation with respect to the value reported since the first ASCA observation (EW~ 200 
eV; Weaver et al. 1997), however when we take into account the larger errors on the early 
ASCA measurements we cannot exclude the line being constant. 

In MCG-5-23-16 this lack of variability of the iron emission line is not so striking as in MCG- 
6-30-15 (Miniutti et al. 2007; Vaughan & Fabian 2004), due to the low level of variability of 
the intrinsic continuum (30-40% compared to a factor of 2-3 in case of MCG-6-30-15). 
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5.3. The blue-shifted absorption line: a possible high velocity, ionized variable 

outflow 

Perhaps the most interesting result of this long XMM-Newton observation has been the 
discovery of a possible variable absorption line from ionized iron. The feature appears to 
be transient with a time scale of about 20 ksec and it is detected at an observed energy 
of about 7.66 keV (corresponding to a rest-frame energy of 7.72 keV). As shown the most 
plausible association of this feature is with K-shell absorption from H-like iron, which is 
blue-shifted by ~ 0.1c. Indeed modeling this absorption feature with XSTAR (Bautista & 
Kallman 2001) requires a column density of about 8 x 10 22 cm -2 and a high ionization param- 
eter (log£ = 3.7 ±0.3 erg cm s -1 ) which implies that the absorption is due to a blueshifted 1 
— * 2 transition of FeXXVl (E — 6.97 keV). The velocity of the absorbing material is found to 
be v = (0.09±0.01)c. In the last few years, red- and blue-shifted absorption lines associated 
with the presence of highly ionized gas in- and/or out-flowing at relativistic velocities have 
been reported for several AGN (E1821+643 Yaqoob & Serlemitsos 2005; Mkn 509 Dadina et 
al. 2005; NGC3516 Nandra et al. 1999; Turner et al. 2005, Mrk 335,Longinotti et al. 2007). 
These absorption lines are found both in Seyfert galaxies (NGC3783, Reeves et al. 2004; 
IC4392a Markowitz et al. 2006; NGC1365 Risaliti et al. 2005; Ark 564, Papadakis et al. 
2007) as well as in quasars (PG 1211+143 Pounds et al. 2003; PDS 456 Reeves et al. 2003) 
and BAL QSOs (APM 08279+5255 Chartas et al. 2002; PG 1115+080, Chartas et al. 2003). 
These absorption systems can also be variable on different time scales, in their ionization 
state and column density, with the most extreme cases being NGC 1365 (Risaliti et al. 2005) 
and Mkn 509 (Dadina et al. 2005). From the analysis of the different intervals in which we 
split the observation we can infer that it is unlikely that the variability of the absorber in 
MCG-5-23-16 is due to a change in the ionization state of the outflowing material, otherwise 
we would detect signatures of this absorber in all the spectral slices. A more likely scenario 
is a change in column density of this absorber. 

The present data suggest we are seeing a transient absorber, which could be associated with 
a possible cloud which sporadically obscures the central source. This “cloud” could be the 
signature of a clumpy absorber located close to the central X-ray source or of matter ejected 
sporadically. Different models have been proposed to explain the powerful outflows detected 
with the recent XMM-Newton and Chandra observations; in particular transient red and 
blue shifted absorption lines are predicted in several theoretical models of failed disk winds 
(Proga et al. 2000; King & Pounds 2003) or an aborted jet (Ghisellini et al. 2004). The pic- 
ture emerging is that these features can provide a direct probe of the dynamics, kinematics 
of the innermost central regions of AGNs. For the absorber detected in MCG-5-23-16 the 
data suggest that this feature appears when the source reaches a relative maximum in the 
intrinsic 2-10 keV flux. However monitoring the spectrum on longer timescales would be 
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required to determine whether there is a statistically firm correlation between the absorber 
and the source brightness or if any duty cicle is present. Thus both a scenario where a 
clumpy absorber or a variable or failed outflow or jet are at present possible. 

Before deriving an estimate of the location, mass and energetics associated with the absorber 
we performed a consistency check for the Nh value measured with the XSTAR model and the 
EW (~ 50 eV) of the absorption line measured with the Gaussian component. Following the 
curve of growth for H-like iron (see Fig. 4 Risaliti et al. 2005) we inferred that the detected 
EW requires a turbulence velocity greater than 500 km/s; a lower turbulence velocity would 
imply a Compton thick absorber and a similar EW of the H-like K(3 line (due to saturation 
of the Ka line), which is not detected (|£W| <15 eV at 90% confidence level). On the 
other hand a turbulence velocity greater than 3000 km/s would produce a broad absorption 
feature which would be resolved even at the EPIC-pn resolution at 8 keV (~ 170 eV). We 
therefore conclude that the observed EW and line width are in broad agreement with the 
XSTAR estimate of a column density of about 10 23 cm -2 , for a turbulence velocity a = 1000 
km/s. 

Using this value for the column density we can now estimate the maximum distance of this 
cloud or blob from the central black hole using the relation between the ionization param- 
eter, the density of the absorber and the illuminating continuum luminosity: L/£ = nR 2 ; 
where L is the intrinsic 2-100 keV X-ray luminosity (5.4 x 10 43 erg/s). Assuming then that 
the thickness of the cloud A R = Nn/n is less than the distance R, we found R < 10 17 cm. 
A lower limit for the distance of the absorber - assuming it is in the form of an outflow - 
can be obtained equating the velocity of the absorbing material to the escape velocity at a 
given radius R from the central black hole; the derived distance is then R £ 100 R g . A con- 
straint on the size of this cloud can be placed assuming 20 ksec as characteristic variability 
timescale, when our line of sight intercepts the aborbing cloud; this gives us A R ~ 6 x 10 13 
cm, which corresponds to ~ 10 R g for a black hole mass of 5 x 10 7 M o (Wandel & Mushotzky 
1986). We can then infer a density of ~ 10 9 cm -3 and assuming a spherical geometry for 
the cloud, a mass of about 10 28 g. These one order of magnitude estimates for the mass 
and velocity correspond to a kinetic energy E kin ~ 5 x 10 46 erg and using the 20 ksec as a 
characteristic timescale to a power of ~ 2.5 x 10 42 erg/s. This value corresponds to ~10% of 
the 2-10 keV X-ray luminosity and is thus is agreement with a radiation driven wind model 
(Proga & Kallman 2004) or with an aborted jet (Ghisellini et al. 2004). 

In conclusion this deep XMM-Newton observation revealed that the soft X-ray emission of 
MCG-5-23-16 can be ascribed to material photoionized by the AGN, likely to be located 
outside the sub pc scale of the absorber and perhaps coincident with the NLR. We confirm 
the presence of a iron Kcc emission line complex composed by a narrow and a broad relativis- 
tic component. The inclination derived from the diskline profile ( i ~ 40°) is in agreement 
with the orientationally based Unification Scheme of AGN (Antonucci 1993) and the X-ray 
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classification of MCG-5-23-16 as Compton thin Seyfert 2 (i.e. intermediate between a type 
1 AGN and a Compton thick type 2) and the optical classification as a Seyfert 1.9. Finally 
we detected a sporadic Fe K absorption feature which could be a signature of a variable high 
velocity outflow. This detection adds one more example to the increasing sample of AGN 
where relativistic outflows have been revealed in the X-ray band. The growing evidence of 
high velocity outflows in AGN indicates that they may play an important role in the ener- 
getics of AGN central engine. 
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Fig. 1. — Ratio of the XMM-Newton pn data and the model when fitting an absorbed and 
redshifted power law over the 2-10 keV band ignoring the iron emission line energy range 
(5.5 — 7.5 keV) . The photon index of the power law has been fixed to 1.8. Two different 
residuals are clearly present: a soft X-ray excess below ~ 2 keV; and the Fe complex at ~ 6.4 
keV 
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Fig. 2. — Ratio of the 0.5— 1.5 keV energy range XMM-Newton pn (black points), MOSl 
(red) and MOS2 (green) data and the model, when the soft X-ray emission is fitted by adding 
to the baseline model, composed by an absorbed and redshifted power law (r = 1.8), a soft 
power-law (r ~ 3) component. An emission line is clearly detected at ~ 0.9 keV (EW~ 60 
eV). 
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Rest Energy (keV) 

Fig. 4. — Ratio between the pn data and the absorbed power-law model (r = 1.65) showing 
the iron line profile. The data clearly show a narrow core at 6.4 keV, red and blue wings 
extending from 5.7 keV to 7 keV and a narrow emission feature at ~7.05 keV, which is due 
to Fe K(3. A sharp drop is also present at ~ 7.1 keV due to presence of reflection. 
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Fig. 5. — Residuals of the data/model of MCG-5-23-16 pn data at the Fe band. Panel a): 
data/model ratio when the underlying continuum is composed of an absorbed power law 
plus a reflection component (R=l.l). Panel b): residuals left when the narrow Fe K a and 
Fe K/3 are added to the model. An excess between 5.8 keV and 7 keV is still present. Panel 
c): residuals when the broad component of the Fe Ko; line is fitted with a diskline model . 
The overall fit is now good and only a weak absorption feature is left at 8 keV. 
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Fig. 6. — Data/model ratio of the Chandra HETG. A strong narrow core at 6.4 keV is 
observed. When modeled with a single Gaussian the line is found to have EW ~ 80 eV and 
a width of a — 32+ eV. 
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Fig. 7. — Contribution to the y 2 in the 5. 5-8. 4 keV for the five pn spectra extracted with a 
time bin of 20 ksec. The Fe Ka line has been parameterized with two Gaussian lines; all the 
parameters of the model except the primary power law normalization have been kept tied 
together. The only strong deviation in the Ay 2 is present in the third spectrum at 7.7 keV. 
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Fig. 8.- — Time resolved spectral analysis. Panel a): Fe 6.4 keV narrow core intensity (in 
unit of photons cm -2 s -1 ) versus the time intervals. The dashed lines correspond to the 90% 
confidence level of the normalization of the narrow core measured in the average spectrum. 
Panel b): same as panel a) for the broad component. Panel c): Absolute intensity of the 
absorption feature (in unit of photons cm -2 s -1 ), the line energy has been fixed to the best 
fit value found for the third interval (7.71 keV). Panel d): MCG-5-23-16 2-10 keV flux (in 
erg cm _2 s _1 ) . Error bars and upper limits are at 90% confidence level. 
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Fig. 9. — Time resolved spectral analysis. Data/model ratio of the pn spectrum of the third 
interval, with the data binned with a constant energy bin of 80 eV. A deep absorption feature 
is visible at 7.7 keV (see text for details). 
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Fig. 11. — Possible high ionization warm absorber models. Panel a) Best fit model: an 
ionization state characterized by log £ = 3.7^oi erg cm s -1 and a column density of ~ 8 x 10 22 
cm ~ 2 with an outflow velocity of 0.09 ±0.01 c. The FeXXVi Lya absorption feature shifted 
to 7.7keV fits well the obseerved deep at this energy. Panel b) The absorption feature at 
7.7 keV can be accounted for assuming a column density of A/fi = 10 23 cm -2 , an ionization 
parameter of log£ = 3.0 erg cm s -1 and no velocity shift. However, the model predicts 
stronger absorption due to the Fe XXIV ls-2p absorption line at 6.67 keV, which is not 
detected. Panel c) A lower ionization Kp absorption line without strong Ka from the 
same species, can be obtained assuming a column density of Nw = 10 23 cm~ 2 , an ionization 
parameter of log £ = 1.5 erg cm s -1 . In this scenario a blueshift of ~ 0.08c would still be 
required to model the absorption line at 7.7keV in the spectrum and the column density 
required to model the EW of the K P absorption feature is Nh >> 10 23 cm -2 , which would 
introduce too much continuum absorption below 6 keV, inconsistent with the pn data. 
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Table 1. Log of the observations and exposure times. 


Mission 

Instrument 

^(total) (k s ) 

T(net) (k s ) 

Tstart 

Tstop 


XMM 

PN 

131.5 

96.2 

08-12-2005T21: 11:28 

10-12-2005T09: 17:08 


XMM 

MOS1 

131.5 

101.6 

08-12-2005T20:42:51 

10-12-2005T09: 16:48 


XMM 

MOS2 

131.5 

102.8 

08-12-2005T20:43:51 

10- 12-2005T09: 16:53 


XMM 

RGS1 

131.6 

97.2 

08-12-2005T20:41:37 

10-12-2005T09: 18:03 


XMM 

RGS2 

131.6 

97.2 

08-12-2005T20:41:42 

10-12-2005T09: 18:03 


Chandra 

ACIS-S HETG 

30 

- 

08-12-2005T17:41:30 

09- 1 2-2005T02 : 33 : 59 


Chandra 

ACIS-S HETG 

20 

- 

09-12-2005T20:52: 1 1 

10-12- 2005T03:00: 10 



N 0 t e . — For XMM- Newton the exposure values reported are total and net exposure time after filtering for high-background 
time intervals. 
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Table 2 . XMM-Newton RGS 1 and RGS 2 : best fit emission lines required fitting with a 

photoionization model. 


Energy 

Flux 

ID 

EW 

AC 

ELab 

(keV) 

( 10 _5 ph cm -2 s -1 ) 

(eV) 

(eV) 


(eV) 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

0.499!°;°°} 

0 . 68 l°;| 

N VIlLyo; 

13.4 

14.3 

0.500 






0.561 (f) 

0.5641“;°“ 

9 oq+2.25 
10 

O vil Hea 

57.2 

22. .5 

0.569 (i) 
0.574 (r) 

0-6531°;“} 

0-631°;! 

O VIII Lya 

20.1 

15.4 

0.654 

n 709 + 0.006 
u. / oz _ 0 005 

0-581°;! 

O VII RRC 

22.7 

6.8 

> 0.739 

0-9001°;“° 

0.91!°;“ 

Neix Hea 

52.3 

10.2 

0.915 (i) 
0.922 (r) 


Note. — The T of the soft power-law tied to the hard power law 
component. The energy of the lines are quoted in the rest frame. Fluxes 
and possible identifications are reported in column 2 and 3. The EW 
are reported in column 4 and they are calculated against the soft power 
law component. In column 5 the improvement of fit is shown using the 
C-statistic; the value for the model with no lines is C=455.7 for 358 PHA 
bins. In column 6 we report the theoretical value for the transitions. 



Table 3. 


Results of the fit for the mean spectrum and the low and high flux states. 



Parameter 

Mean 

High 

Low 

Continuum 

r 

1.82 ±0.01 

1.84 ±0.01 

1.81 ±0.01 


N| 

1.49 ±0.01 

1.50 ±0.02 

1.49 ± 0.02 


Flux 6 

8.16 

8.84 

7.39 


Narrow Gaussian 

E 

6.42 ± 0.01 

6.41 ±0.02 

6.40 ±0.02 


EW 

611? 

48ig 

811 “ 


N c 

5.4±8;1 


6-511:1 

Broad Gaussian 

E 

6.4 f 

6.4 f 

6.4 f 


a 

0.35 ±0.1 

0.35l°;l 6 

0.3712; 22 


EW 

64 

61^23 

71X+32 

'°-32 


N c 

5.9 ±1.5 

5.91& 

6.3llj 

Diskline 

R-in 

48^20 




i 

4iig 




EW 

co+14 

00 -13 




N c 

4.6 ±1.2 




Note. — The line energies are expressed in units of keV, while their a 
and EW in eV. 

a Column density in unit of 10 22 cm -2 
b 2-10 keV flux in units of 10 -11 erg cm -2 s -1 
Normalization of the Fe line in unit of 10 -5 photons cm -2 s -1 
indicates that the parameter has been kept fixed. 
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Table 4. Best fit parameters for the absorpion feauture detected at 7.7 keV for the EPIC 

pn and MOS cameras and SUZAKU. 



(keV) 

(10 5 ph cm 2 s *) 

(eV) 

EPIC-pn“ 

7.9 ±0.1 

-1.9 ±0.6 

33tt 

EPIC-pn 6 

7.7 ±0.1 

-3.2 ±0.9 

52 ±15 

EPIC-MOS“ 

7.4 ±0.2 

-1.1 ±1.0 

16 ±14 

EPIC-MOS 6 

7.4 ±0.2 

-2.2 ± 1.8 

27 ±23 

SUZAKU-XISO- 1-2-3“ 

7.8 ±0.1 

-1.8 ±0.9 

30 ±15 

SUZAKU-XIS0-1-2-3 6,c 

7.65 ±0.4 

-3.5 ±3.0 

52 ±44 


Note. — The parameters of the absorption line are derived adding 
an inverted gaussian component to the best fit model. 

a Parameters refers to mean spectrum 

b Parameters are derived for the 20 ksec time where the line is 
detected. 

c The net exposure of this Suzaku spectum is only ~ 6 ksec 



